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Abstract 

Rapid oscillations in the corona will be discussed from a theoretical 
standpoint, developing some previous work on ducted, fast magnetoacoustic 
waves in an inhomogeneous medium, in the theory, impulsively (e.g. flare) 
generated mhd waves are ducted by regions of low Alfv£n speed (high density) 
such as coronal loops. Wave propagation in such ducts is strongly 
dispersive and closely akin to the behaviour of hove waves in seismology, 
Pekeris waves in oceanography and guided waves in fibre optics. Such flare- 
generated magnetoacoustic waves possess distinctive temporal signatures 
consisting of periodic, quasi-per iodic and decay phases. The quasi- per iodic 
phase possesses the strongest amplitudes and the shortest time scales. Time 
scales are typically of the order of a second for inhomogeneities (coronal 
loop width) of 1000 km and Alfven speeds of 1000 Jans -1 , and pulse duration 
times are of tens of seconds. 

Quasi-per iodic signatures have been observed in radio wavelengths for 
over a decade and more recently by SMM. It is hoped that the theoretical 
ideas outlined above may be successfully related to these observations and 
thus aid the interpretation of oscillatory signatures recorded by SMM. Such 
signatures may also provide a diagnostic of coronal conditions. 

New aspects of the ducted mhd waves, for example their behaviour in 
smoothly varying as opposed to tube- like inhomogeneities, are currently 
under investigation. The theory is not restricted to loops but applies 
equally to open field regions. 
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1. Introduction 


Since the end of the 1960's many observations of coronal oscillations 
have been reported, the bulk of the evidence being provided by radio wave 
and hard X-ray data (Abrami, 1970, 1972 j Rosenberg, 1970) McLean et al . 

1971; Gotwols, 1972; Kai and Takayanagi, 1973; McLean and Sheridan, 1973; 
Kaufmann et al . 1977; Pick and Trottet, 1978; Tapping, 1978; Dennis et al . 
1981; Trottet et al . 1981; Kane et al . 1983, Kiplinger et al . 1983 a,b; 
Takakura et al . 1983), though coronagraph observations have also been 
reported (Koutchmy et al . 1983; Pasachoff and Landman, 1984). Besides the 
long period (one minute and greater) radio pulsations reported by, for 
example, Trottet et al (1979) there are numerous observations of short 
period oscillations with periodicities of around one second. An interesting 
feature of these short period observations is the form of the observed 
signature, the temporal variation. Short bursts of oscillatory behaviour 
are seen to occur, the wave packet or wave train behaviour being described 
variously as periodic, quasi-per iodic or as a sequence of spikes. Kiplinger 
(1984, private communication) has pointed out that of the various criteria 
that any successful model of these oscillations must meet, one is certainly 
that the process is not am enduring one. Theory must be able to explain the 
termination of the pulses. Kiplinger also stresses two other observed 
features of the hard X-ray (and also microwave and y-ray) pulse train 
phemonenas the source appears to be more or less stationary and compact, 
and the oscillations are very irregular. A model should ascribe a time 
scale of the variations rather than a well-defined period i.e. one phase of 
the motion is quasi-periodic . 

How are these rapid oscillations generated? It is suspected that the 
source is impulsive, perhaps a solar flare, as in the examples given by 
McLean and Sheridan (1973) and Kiplinger et al (1983 a,b). However, Tapping 
(1978) and Gaizauskas and Tapping (1980) have reported bursts of 
oscillations with which there appears to be no associated flare activity. 


2. Model of a Coronal Structure 

The open and closed 'loops' of the magnetically dominated corona vary 
in size and are often twisted, hot and dense. Structures may have cool 
cores and it is thought that brightness variations over them may imply 
density variations . Not all these observed properties can be incorporated 
into a model Which is simple enough to exhibit the wave-like phenomena 
associated with the corona. Here the effects of gravity (not so important 
for a region of scale >, 3 x lo*km above the photosphere), curvature (more 
important for closed loops ) and twist ( important in knowing if an 
equilibrium exists) will be neglected. Since transverse dimensions are 
usually much smaller them longitudinal ones, the structure will be modelled 
by a straight cylinder and will be regarded as a density inhomogeneity. 
This may take the form of a uniform magnetic cylinder or alternatively the 
density may vary across the inhomogeneity but in either case we consider 
this to be an intrusion into an otherwise uniform magnetic environment. 

It has been shown previously (Roberts, Edwin and Benz, 1984) that a 
uniform, high density coronal cylinder of radius a, its axis lying along the 
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z-direction, acts as a wave guide for the free nodes of oscillation 
(e i (* z+w t>, u, k real) of the structure. A high density (v^ > v A , where 
v A , v Ae 8X6 the Alfven speeds in the cylinder and its environment 
respectively) coronal loop can oscillate freely (that is without radiating 
energy to infinity) in the fast magnetoacoustic mode. Motions in the 
environment of the coronal loop arise simply in response to those within the 
loop i the loop acts as a wave duct trapping the fast mihd oscillations. Which 
penetrate only a small distance into the environment. 

The dispersion relation for these ducted, free modes of oscillation is 
(Meerson et al . 1978; wentzel 1979; Wilson 1980; Spruit 1982; Edwin and 
Roberts 1983; Roberts et al . 1984) 
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and p Q and p e cure the equilibrium densities inside and outside the cylinder, 
respectively. The sound and Alfven speeds inside the cylinder are c Q and 
v A , outside c e and v^. Also, 
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Kjj, J n are Bessel functions (modified and first kind, respectively), n - 0 
give the symmetric (sausage) modes (those in which the axis of the cylinder 
is undisturbed by the vibration ) and n — 1 the kink or asymmetric modes . 

Figure 1 illustrates the solutions of Equation (1) for a dense coronal 
loop: two sets of waves feature which, because the sound and Alfven speeds 
differ by an order of magnitude, have widely separated frequencies. One 
class of oscillations has slow, acoustic time scales, for both kink and 
sausage inodes, and the waves sure only mildly dispersive. In a low 0 plasma 
the tube speed c^, is close to the sound speed c Q so the frequencies of both 
these slow modes are given (to a good approximation, provided ka is not too 
large ) by 

w - kQp » kc Q . (2) 

Of the fast waves (with phase-speeds satisfying v A < u/k < v^ ) only 
the principal kink mode exists for all wave numbers k (i.e. does not have a 
cut-off). It has a typical frequency given by 

“ - k< V (3) 


where 
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Figure 1 The phase speed u/k (u, k real) as a function of ka (k > 0) 

where a is the radius of a dense coronal structure, illustrated 
here for v^ — 5c Q , c e " 0.5c o and v A = 2c Q . The mildly 
dispersive, slow band of oscillations has infinitely many modes, 

on acoustic time scales, for both sausage ( ) and kink ( ) 

modes. Here only two cure shown. Except for the principal kink 
mode, which exists for all wave numbers k, the fast waves have 
cut-offs. Hatching denotes regions where there are no free modes 
(After Edwin & Roberts, 1983). 


In a coronal loop of length L, with its footpoints anchored in the 
dense chromosphere/photosphere both the slow and fast waves may occur as 
standing modes. Then Equation (2) and (3) would give rise to two time 
scales T g - 2L/c 0 and r f - 21/c^, corresponding to the slow and fast waves 
respectively. For -typical coronal parameters, for example a sound speed of 
200 km s -i , an Alfven speed of ~ 2000 km s -1 and a loop length of 2 x io 4 - 
10* km, t 8 and r f are estimated to be of order 200- 1000s and 15-708 
respectively. 


3. Impulsively Generated Modes 

A solar flare, considered as either a single or multiple impulsive 
source, will initiate propagating rather than standing waves. 
Alternatively, there may be other less energetic sources generating the 
waves. Which may arise either in a closed coronal loop, if the motions have 
insufficient time to reflect from the far ends of the loop, or in open field 
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regions . 

To simplify the discussion and interpretation of these propagating 
waves it will be assumed that the plasma is cold. This is reasonable for 
the low 0 corona and involves no loss of essential physics an far as the 
fast inodes axe concerned. Mathematically it means that c Q - c e - o in 
Figure 1 and that just the fast band of waves remains. For the fast sausage 
waves in particular the dispersion relation (1) with c e = c Q = 0 reduces to 


n K i.( m e a > _ J x< n o a > 

° W 1 * ^ J o< n o a > 
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where 
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For a cylindrical structure of large radius , ( 4 ) reduces to 

tan(n 0 a - ir/4) - - (5) 

Which, except for a phase shift of it/4, is identical to the dispersion 
relation given by Pekeris (1948) for the propagation of explosive sound in 
shallow water. 

Thus Eqn (4) describes the highly dispersive, fast, symmetric (cold 
plasma) waves of Figure 1. If such waves are generated impulsively, then 
the resulting disturbance 4> may be represented as a Fourier integral over 
all frequencies u and wave numbers k. In general, a wave packet results. 
For example, if o> is related to k from (4) by u = ±W( k ) , then W'(k) describes 
the group velocity of the waves and 4> has the form (Pekeris, 1948; Whitham, 
1974 )i 


A(z,t)e i 'W z,t ), W'(k) * 0, 

1 e 1 ** 55 '*), W'(k) - o, 
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(6a) 
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where 


A(z,t) " , e-< ffi /4)sgnw'(k) / ^ * (z , t) . zk _ w(lt)t . 

y t iw ( x ) i 

Thus Eqn (6a) describes an oscillatory wave train with vj» giving the 
variation between local maxima and minima. However, since A, k and w are 
not constant the wave train is not uniform. The distance and time between 
successive maxima, and the amplitude, are not constant i.e. there is a quasi- 
periodic behaviour. It is obvious from Eqns (6) that extrema, W'(k) = 0, of 
the group velocity are important since then expression (6b), the so called 
Airy phase, describes the disturbance rather than (6a). 

The group velocity curve for the m-k dispersion relation (Eqn (4)) for 
a uniform cylinder of density inhomogeneity is found to possess a minimum. 
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So too does the group velocity curve for a ji-power density intrusion into an 
otherwise uniform medium. If the density variation across the loop has the 
form 


Po< r > - 


Po<°> 


■ ffll 1 '£]]• 


0 < r < a. 


r > a. 


(7) 


then the u-k dispersion relation becomes (Edwin, 1984) 


v K a (v) 

K A (v) 


K 0 (v) K 2 (v) 


Kf(v) 


tl+2 


u J Q (u) 

J X (U) 


( 8 ) 


where u 2 * — ^ — n 2 a 2 and v 2 - m£a 2 . Letting ii - ® in the relation (8) 

H + 2 ° e 

recovers the uniform density dispersion relation Eqn (4). 

Typical group velocity plots for the lowest order sausage modes of the 
uniform (n - oo) and parabolic (n - 2) density variation cylinders axe shewn 
in Figure 2. 



Figure 2 The dimensionless group velocity c_/v Ae as a function of 
dimensionless frequency iaa/v Ae of the lowest order fast 

sausage mode of Figure 1, shown for uniform ( ) and parabolic 

( ) density profiles across a coronal structure. Here p Q (0), 

the density on the axis of the structure is 6 times that in the 
exterior of the loop, p e . (After Edwin, 1984). 
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Since each curve has a minimum there is a point where w'(k) * 0. The 
disturbance resulting from an impulsively generated fast symmetric mode may 
thus be described as follows. At time t - 0 a pulse is generated at a point 
source z = 0, say, which is composed of all frequencies. The disturbance 
observed at a (large) distance z - h from the source, on the basis of one of 
the group velocity curves in Figure 2, initially has a frequency given by 




*c v Ae' 


( 9 ) 


Where the wave number k c at cut-off is found from (4) or (8) with w/k * Vj^. 
Specifically, 


f.j<°> 

— — , S— 1,2,3,... 

(P o (0)/p e - l )*a 


( 10 ) 


k 

where f - 1 and p o (0) - p Q for the uniform cylinder, and f = (1 + 2 /p) for 
the p - power density variation! (-2.40, 5.52, ...) are the zeros of the 

Bessel function J Q . nils disturbance has taken a travel time of t - h/v Ae 
from the source and initiates the periodic phase. As time progresses the 
frequency and amplitude of this nearly sinusoidal wave (Eqn (6a)) grow 
slowly until a time h/v A when high frequency information arrives from the 
source. So a new train of high frequency waves, due to the right-hand branch 
of the group velocity curve, is superimposed on those arriving from the 
left-hand branch. The result is a strong increase in amplitude and the 
oscillation becomes quasi-perlodic (see (6a) where the amplitude varies 
inversely as the square root of the slope of the group velocity curve). 
During this quasi-perlodic phase the frequencies of the waves in the 
right-hand and left-hand branches of the curve continue to decrease/ increase 
respectively (see Figure 2) until at a time t = h/Cg in , where Cg in is the 
minimum group velocity, they coincide. The disturbance then consists of the 
single frequency u min (Figure 2) and, though the disturbance continues to 
oscillate with this constant frequency, the amplitude declines rapidly 
according to (6b) (W(k) - <J nin ) in the Airy or decay phase. The behaviour at 
a location z - h from the impulsive source is well known (Pekeris, 1948; 
Ewing et al . 1957; Brekhovskikh, I960; Tolstoy, 1973; Kennett, 1983) and the 
sequence of phases just described can be illustrated as in Figure 3 . 

In order to draw comparisons between observations and theory the 
particular frequencies and time scales associated with Figure 3 must be 
identified. The frequency representative of the periodic phase is given by 
(9) and (10), and setting s - 1 for the lowest order sausage mode. In terms 
of period this means that the largest periodicity of the disturbance is given 
by 


2tra 




'Ae 
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For a very dense inhomogeneity (p Q >> p e ), Eqn (11) reduces to 
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periodic phase 


quasi-periodic decay 
phase phase 



Figure 3 & sketch of the evolution of the lowest order fast sausage 

node of Figure l in the low 0 extreme (c Q , c e << v A , v^ ) 
showing the various phases in the disturbance as recorded 
at an observation level z - h away from an impulsive source 
at z * o. 



( 12 ) 


The other important frequency is that representative of the quasi- 
periodic phase, though this is a variable quantity because of the very 
quasi-periodicity. u ia ^ n , which characterizes the frequency at the end of 
the quasi-periodic phase is larger than « c and so has an associated period 
* m * n Which is smaller than t c . Thus rapid oscillations occur during the 
quasi-periodic phase . 


4. Relation to observations 

The mechanism by which fast magnetoacous t ic waves ducted in a density 
enhancement of the corona may be manifest as radio wave observations was 
discussed in Roberts et al (1984). Here we merely summarise the properties 
of such mhd waves when they are inpuls ively generated, perhaps by a solar 
flare, and show that several of the desirable features of observed rapid 
oscillations can be explained by our model. 

When the lowest order, fast, symmetric mode is generated impulsively it 
exhibits both periodic and quasi-periodic signatures as shown in Figure 3, 
the vast majority of the wave power appearing as a well-defined packet. The 
time signature thus resembles the pulsations recorded in the radio wave (and 
other) data (see e.g. Pick and Trottet, 1978; Tapping, 1978). The time 
scales associated with such a disturbance are indeed rapid and of the order 
of 1 second, as observations require. For example , taking a particle 
density of ~ 10®cm -3 , a magnetic field intensity of 40G and a loop radius of 
- 10®km, Eqn (5) shows that the periodicity at cut-off for a uniform loop is 
0.94s. Again it must be remembered that this is the largest time scale of 
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the motion, the periodicity within and at the end of the quasi-per iodic 
phase being much smaller. Also coronal loops having a gradual density 
variation with radius rather than forming a step-like density inhomogeneity, 
would according to the ji-power model have even larger values of and a 
correspondingly reduced time scale of oscillation in the quasi-periodic 
phase. 

Apart from periodicities, there axe time scales associated with the 
durations of the periodic and quasi-periodic phases. observational 
determination of the onset of the periodic phase is expected to be difficult 
because the lower amplitude oscillations may be lost in background noise, 
but the quasi-periodic phase is theoretically distinctive (Figure 3) and so 
comparisons with estimates of its duration, T dur , given by 



and observed signatures of similar form should be possible. It was 
suggested also in Roberts et al . (1984) that the width of the irihomogeneity 
and the associated Alfven speeds may be obtained by plotting a/h (in units 
of where T onset ” h / v A> against T dur (in units of r onset ), so 
that the theory may be used to determine physical conditions in the corona 
using radio pulsation data. 

The picture presented by Figure 3 depends largely on the comparative 
sizes of v A , v^ and Cg in . If the density of the ihhomogeneity is not very 
different from that of its surroundings, so that v A is close to v^ and the 
structuring is weak, then the quasi-periodic phase is of fairly short 
duration because cJJ* 11 - v A (see Figure 2). Consequently, with the periodic 
phase lost in noise, the observed pulse train will start abruptly. Also it 
will decay more gradually because the Airy phase is correspondingly 
extended. (The curvature of the group velocity curve at the minimum is not 
so pronounced for v A - v^ so |W"(k)| in (6b) is smaller.) This theoretical 
possibility may correspond to the event observed by McLean and Sheridan 
(1973). Density variations across the loop would produce a lower c® 1 ” and 
longer duration times (see Figure 2 and Eqn (13)). 


5. Concludi ng BpmarVH 

A mechanism has been offered to try to explain some of the quasi- 
periodic, rapid, oscillatory behaviour observed in the solar corona. It has 
been shown that the explanation does not depend on the coronal loop having 
uniform density. Variations in density across the loops would modify the 
duration of the quasi-periodic phase and the time scale of the periodic 
phase, but observations and results from radio wave, hard X-ray and other 
data could still be accounted for. 

Most of the foregoing analysis (the group velocity curves, etc) has 
been carried out for both sausage and kink modes. (For the asymmetric kink 
mode in a wide cylinder ka >> 1, the phase- shifted love's equation 

tan(n Q a - ir/4) - «fe/n 0 . 
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describing seismo logical waves ducted in the Earth's crust, is obtained 
instead of Eqn ( 5 ) . ) All inodes would be generated by an impulsive source so 
a combination of the above features is to be expected. Moreover 
disturbances could be generated by a series of Impulses at various times so 
that a quas i-per iodic phase resulting from one disturbance may interact with 
the periodic phase of a subsequent impulse to give a much more complicated 
picture than that presented in Figure 3. Interactions between impulses and 
the various combinations of modes could therefore produce a feu: more 
distorted signature than that described here for a single mode and a single 
impulse. 

Whilst it is recognised that the real situation in a coronal loop is 
unlikely to be that of a single, symmetric mode in a dissipationless, linear 
system arising from a single impulse, the theory nevertheless offers an idea 
which cam bo tested against further observations. The marriage of 
observation and theory offers the possibility of a fruitful diagnostic of In 
situ coronal conditions. 
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